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A new parametric equation of state and quark stars 
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Abstract It is still a matter of debate to understand the equation of state of cold supra-nuclear matter in 
compact stars because of unknown non-perturbative strong interaction between quarks. Nevertheless, it is 
speculated from an astrophysical view point that quark clusters could form in cold quark matter due to strong 
coupling at realistic baryon densities. Although it is hard to calculate this conjectured matter from first princi- 
ples, one can expect the inter-cluster interaction to share some general features to nucleon-nucleon interaction 
successfully depicted by various models. We adopt a two-Gaussian component soft-core potential with these 
general features and show that quark clusters can form stable simple cubic crystal structure if we assume the 
wave function of quark clusters have a Gaussian form. With this parameterizing, Tolman-Oppenheimer-VolkofF 
equation is then solved with reasonable constrained parameter space to give mass-radius relation of crystalline 
solid quark star. With baryon densities truncated at 2no at surface and range of interaction fixed at 2fm we 
can reproduce similar mass-radius relation to that obtained with bag model equations of state. The maximum 
mass ranges from ~ 0.5Mq to > 3Mq. Observed maximum pulsar mass (> 2Mq) is then used to constrain 
parameters of this simple interaction potential. 

Key words quark star, solid quark matter, mass radius relation, massive pulsar 
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1 Introduction 



An understanding of cold quark matter is both 
one of the most challenging problem in particle 
physics and a prerequisite to understand the true na- 
ture of pulsars and pulsar-like objects. However, due 
to both non-perturbative nature of strong interaction 
at low energy and complexity presented by the quan- 
tum many-body problem, it is almost impossible to 
understand such state theoretically from first princi- 
ples. 

Over the decades various approaches to bypass 
these difficulties has been developed, both perturba- 
tively such as that of color super-conductivity Q and 
non-perturbatively such as lattice QCD and QCD- 
based effective models. On the other hand, it has 
been conjectured that quark matter could be in a 
solid state at extreme low temperature present in pul- 
sar interior. This possibility could combine naturally 
with several previous works suggesting the possibil- 
ity that deconfined quark matter might contain quark 
clusters of 3 AT valence quarks into a reasonable 



conjecture that quark clusters could form crystal lat- 
tices. 

Because of the difficulty to obtain detail of the in- 
teraction between quark clusters and therefore equa- 
tion of state of cold quark matter at a few times nu- 
clear densities, it is interesting to apply simple phe- 
nomenological models. If we can use astronomical 
observations to constrain parameters in such models, 
we will be able to gain some insight on properties of 
low-energy QCD or rule out such form of cold quark 
matter within pulsars. 

In several previous works different models 

have been tried to investigate the possible equation of 
state of solid quark matter and have provided possi- 
bility to explain stiffness in equation of state required 
by observed massive pulsars [8| . In Q the Lennard- 
Jones potential which was introduced to model in- 
teraction between inert gas molecules |9| is used as 
potential between two quark clusters. The Lennard- 
Jones potential shares some basic properties with nu- 
clear forces such as short-range repulsion and longer 
range attraction (c.f. [13]). In this article we adopt a 
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more realistic parametrization that shares great simi- 
larity with various models depicting hyperon-hyperon 
potential. It has been shown that the interaction be- 
tween two H-dibaryons - cluster of 6 valence quarks 
may also share this general feature 
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This article is arranged as follows. The model 
of inter-cluster potential is presented in Section 2. 
Parameter space used for calculation is discussed in 
Section 3. Section 4 shows result of calculation. Con- 
clusion and some discussions are presented in section 
5. 

2 Inter-cluster potential 

As in Q we consider quark clusters with 3N va- 
lence quarks with N — 1,6 and a simple cubic lattice 
structure is adopted for simplicity. In the context 
of strange quark matter, these are particles with the 
same valence quark composition as hyperons and the 
hypothetical 'quark-alpha' 12j. By extension of the 



Lennard- Jones potential, we adopt the following sim- 
ple parametrization for effective interaction between 
two quark clusters localized on crystal lattice sites, 



v(r) = V 1 e 



te) 2 



-V 2 e 



(1) 



where possible spin-dependent interactions are omit- 
ted for simplicity. As is mentioned above, with the 
condition V 1 >V 2 , r 1 <r 2 , this potential qualitatively 
reproduces the general feature of various successful 
phenomenological potentials of nuclear interactions: 
soft-core repulsion at short range and attraction at 
longer range. 

It turned out that maximum mass and mass- 



radius curve is not sensitively depend on the value of 
r 2 . Therefore it is reasonable to fix r 2 at 2fm which 
is a typical range of nuclear force. 

For simplicity we assume that cluster center-of- 
mass has a Gaussian wave packet with width w as 
wave function 



1 



7r 3/4 w 3/2 



(2) 



In 0] it is assumed that potential well created by 
surrounding clusters with Lennard- Jones interaction 
is deep enough to trap quark clusters in the potential 
well. To ensure that soft-core potential with given 
parameter can also achieve this we adopt variational 
method, i.e. to determine the value of w by min- 
imizing total energy of a single cluster which is a 
sum of kinetic energy of the wave packet and po- 
tential energy contributed by surrounding cluster lat- 
tices. The result shows that with the range of of pa- 
rameters considered in this work, the width of wave 
packet is rather small compared to inter-cluster dis- 
tance and hence it makes sense to speak of this sys- 
tem of clusters as quark clusters trapped in periodic 
lattice. With this small width, the overlap between 
adjacent wave packets is negligible. Thus it is rea- 
sonable to omit the difference between fermionic and 
bosonic quark clusters. 

To calculate total contribution to single cluster 
potential energy, a sum is taken over a cube of 21 3 lat- 
tices centered around the quark cluster under consid- 
eration. The size of the cube is enough since cluster 
number density in this work will not exceed ~ 10n . 
Hence, twice the total contribution to potential en- 
ergy for a single cluster is 



rv>=( E E E 

\fc 1 = - 10 fc 2 = - 10 fc 3 = - 10/ 



,1/3 



(3) 



where the prime means that the sum omits ki — k 2 — 
k 3 = 0, and 

v(w;r)= J d 3 r'r , w (r')v(r')^ w (r'), (4) 

is the expectation value of potential energy between 
two clusters. 

The total energy density is then 

3 9, 



e = -V(n, w) + nm + + - 

2 Amw 2 8 



(6n 2 y/ 3 hvn 4/3 , (5) 



where the third term comes from contribution of ki- 
netic energy and w is treated as a function of number 
density n. Same as in [7j the fourth term comes from 



zero-point energy of phonon in Debye's approxima- 
tion, with 



— -- — + — ) 

v 3 3 \ vf, V s , I ' 



(6) 



and i>|| , v ± stands for sound velocity of longitude and 
transverse modes respectively. Following the argu- 
ment of [tJ we take v = c/3 because its value have 
only very small influence on the final result. 

Ignoring effect of temperature the pressure can be 
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derived as 

2 d / e \ n 2 <9V n 2 dV dw 3 du; 
dn\nJ 2 dn 2 chu 2mw 3 dn 

(7) 

With energy density and pressure one can establish 



dn 
dr 
dM 
dr 
M(0) 

where n c is central number density of quark cluster. 
Because quark matter is usually expected to be self- 
bound at zero pressure (for instance under bag model 
equation of state), it is reasonable to simply adopt a 
truncation baryon number density n sur f at the sur- 
face. In this work we adopt n sur f = 2n where n is 
baryon number density of normal nuclear matter. 

3 Parameter space 

As is stated above, r 2 is fixed to 2fm. Therefore, 
we have 4 free parameters: quark cluster mass m; 
height of the two Gaussians in the potential Vl,V 2 ; 
range of repulsive core r { . It is appropriate to ex- 
pect that the depth of attractive part of the potential 
might be of the same order of magnitude as typi- 
cal potential between two nucleons in nuclear mat- 
ter. Therefore we fix V 2 at 50MeV and lOOMeV 
respectively. On the other hand, we fix the value 
of m to lGeV and 6GeV for 3-quark clusters and 
the aforementioned hypothetical 'quark-alpha'. Here 
mass of 3-quark cluster is take from mass of A hy- 
peron (1115MeV). The 18-quark cluster - the 'quark- 
alpha' is assumed to contain 6 quarks of each flavor 
which may have a mass less than 6m A but we omit 
this possibly small difference for simplicity. Thus all 
that are left is V\ and r 1 . We adopt a condition, 

V ir l>V 2 rl, (11) 

which ensures that potential energy is always positive 
(i.e., repulsive) when density is very high. 

4 Results and conclusions 

With the above settings of parameter space, we 
drew 4 contour plots of maximum mass calculated for 
4 different sets of (m,V 2 ) which are shown in FigQ] 



equation of state and solve Tolman-Oppenheimer- 
Volkoff (TOV) equation with varying central density 
to get the mass-radius relation. In practice, it is more 
convenient to skip (numerical) determination of equa- 
tion of state and write the TOV equation in terms of 
n(r) and M(r), 

(8) 

(9) 
(10) 



tings with maximum mass exceeding 1.9Af Q are also 
shown (Fig J5J). From FigfT] we can see that for our 
simple soft-core parametrization, maximum mass can 
range roughly from below 1M to about ~ 3M Q for 
solid quark stars with 3-quark cluster forming crys- 
tal lattice and from below 0.5M© to about 2.1M 
for solid quark stars made up of 'quark-alpha' parti- 
cles. On the other hand typical mass-radius relations 
shown in Figj2] are very similar to those calculated 
within bag model EoS (shown as gray curves in Fig 12 
adopted from model SSI, SS2 in [13J) and M-R curves 



= - J3 -^(P(„) +e( „))( W+ 4, t r»P(n))(^) 

= 4irr 2 e(n) 
= 0, n(0) = n c 
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Fig. 1. Maximum mass in unit of solar mass 
of solid quark star. Upper left: m — lGeV, 
V 2 = lOOMeV; Upper right: m = lGeV, V 2 = 
50MeV; Lower left: m = 6GeV, V 2 = lOOMeV; 
Lower right: m — 6GeV, V 2 = 50MeV. Bound- 
ary of contour plots are moved a little upper 
right than Vix? = V 2 r 2 to avoid parameters 
that lead to large error in numerical calcula- 
tions 



Typical mass-radius relation curves of these set- 
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calculated in [7(. This shows that at least for some 
region of parameter space our simple parametrization 
can also produce heavy maximum mass supported by 
observed value of 2M Q . Inversely, with current obser- 
vation we can already restrict parameters in this very 
simple model with only 4 parameters. For instance, to 
get maximum mass larger than 1.9M for m = lGeV 
and V x < 6GeV, we have to restrict r x to below 0.75fm 
when V 2 = lOOMeV and restrict r 1 < 0.6fm and for 
m = 6GcV it requires r 1 < lfm and Vi > 8GeV. 




2 4 6 8 10 



Radius/km 

Fig. 2. Mass-radius relation curves for typi- 
cal parameters giving large maximum masses. 
Solid: m = lGeV, Vi = 3.2GeV, V 2 = lOOMeV, 
n = 0.68fm; Dashed: m = lGeV, Vi = 5GeV, 
V 2 = 50MeV, ri = 0.54fm, Dotted: m = 6GeV, 
Vi = 9GeV, V 2 = lOOMeV, n = 0.95fm; Gray 
solid and dashed: Bag model EoS SSI and SS2 
adopted from [13|] 

5 Discussion 

In cold quark matter at baryon number densi- 
ties realistic for compact stars, the interaction be- 
tween quarks could be strong enough that instead 
of condensating in momentum space to form color- 
superconductive phase it is possible that dressed 
quarks undergo condensation in position space to 
form quark clusters. As is stated in [7| if the po- 
tential well formed by neighboring clusters are deep 
enough to trap quark cluster, cold quark matter could 
form crystal solid in low temperature. 



In this work we discussed simple-cubic lattice 
structure formed by 3 and 18-quark clusters using 
a simple two Gaussian component parameterization 
of soft-core potential to simulate the interaction be- 
tween quark clusters. This parametrization shares 
the basic properties of nucleon-nucleon interaction 
mediated by meson exchange - short range repulsion, 
medium and long range attraction and a finite range. 
These properties are also shared by Lennard-Jones 
potential adopted in Q. However unlike Lennard- 
Jones potential with r~ 12 pole at origin, soft-core po- 
tential adopted in this work can be treated by non- 
relativistic quantum mechanics. By minimization to- 
tal energy we found that at realistic densities this 
soft-core potential can lead to a stable lattice struc- 
ture. It is entirely possible that other unit cell struc- 
ture (e.g. body-centered cubic) is more stable, but 
we expect the difference to be quantitative instead of 
qualitative. 

It is also worth mentioning that despite great simi- 
larity between mass-radius relation obtained with bag 
model equation of state and those calculated here, 
the underlying picture is drastically different. In bag 
model, quark star is degenerate Fermi gas of free 
quarks sustained by vacuum energy and associated 
negative pressure. In this work and 0] pressure is 
mainly provided by repulsive core of inter-cluster po- 
tential instead of mere degenerate pressure. 

Although maximum mass cannot be easily tuned 
to ~ 6M Q as in J] due to soft-core nature of the 
interaction, our parametrization can still provide a 
stable lattice crystal structure with maximum mass 
exceeding 2M which is in accordance with observed 
maximum massQ. Inversely, the observed maximum 
mass can be used to put constraints on parameters 
of this simple model which will possibly give some 
insights into the form of interaction between quark 
clusters if such phase exists. 

We thank pulsar group at Peking University for 
useful discussions. 
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